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Abstract 



We analyze previously proposed order parameters for the confine- 
ment - deconfinement transition in lattice SU(2) Yang-Mills theory, 
defined as vacuum expectation value (v.e.v.) of monopole fields in 
abelian projection gauges. We show that they exhibit some incon- 
sistency in the treatment of small scales, due to a violation of Dirac 
quantization condition for fluxes. 

We propose a new order parameter avoiding this inconsistency. It 
can be interpreted as v.e.v. of the field of a regular monopole in 
any abelian projection gauge, but it is independent of the choice of 
the abelian projection. Furthermore, being constructed in terms of 
surfaces of center vortices, it has also a natural interpretation in the 
approach of center dominance. 



1. Introduction 

There now appears to be a general consensus about the idea that color con- 
finement in Yang Mills theories is caused by the condensation of some topological 
defects. However it is still debated whether the important defects are center vor- 
tices or magnetic monopoles. 
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The first proposal was brought forward, in the '70, by 't Hooft, Polyakov, Mack, 
Nielsen and Olesen [1] and others and has received renewed interest in recent times, 
stimulated by numerical simulations [2]; presently it is usually named "center 
dominance" . 

The basic idea is drawn from an analogy with ferromagnets and may be roughly 
phrased as follows. The expectation value of a temporal Wilson loop in four 
dimensions can be seen, by fixing the temporal gauge, as a product of two-point 
spin correlation functions of three dimensional non linear a models coupled to 
gauge fields, one for each time involved in the Wilson loop. If the non-linear a— 
models are disordered, their spin correlation functions have an exponential decay. 
This implies an area law for the Wilson loop of the four-dimensional gauge theory 
and hence, according to Wilson's criterion, confinement. (Actually, a rigorous 
proof of a mathematically precise version of this idea has been established in [3]). 
The defects responsible for disorder in the non linear a models are vortices. For 
SU(N) theories they are Z^r center vortices. 

The alternative with monopoles as the relevant defects [4] has been put on a 
concrete basis by 't Hooft [5], who suggested to make explicit the monopoles in 
SU(N) Yang-Mills theories, by performing a gauge fixing that leaves a maximal 
Cartan torus, {/(l)-^ -1 , unbroken. These gauges are called abelian projection 
gauges. The resulting abelian gauge theories can be rewritten in terms of "pho- 
tons", gauge fields of a theory with gauge group the decompactification of the 
residual gauge group, and monopoles. The monopoles corresponding to this de- 
compactification are those expected to condense in the confinement phase. 
To define in SU(2) Yang-Mills theory, to which this paper is addressed explicitly, 
an abelian projection gauge, 't Hooft suggested to construct some scalar field X 
taking values in su(2), as a function of the gauge connection and transforming in 
the adjoint representation of the gauge group, SU(2). To perform the gauge fixing 
one imposes the constraint that X is diagonal. 

The diagonal component of the 5 , L r (2)-connection in this gauge plays the role 
of the "photon field" , the off-diagonal components are charged with respect to the 
residual gauge group U(l). 

The points in space-time where two eigenvalues of the matrix X coincide are the 
positions of the monopoles in this gauge. Confinement is then believed to emerge 
as a consequence of monopole condensation in the form of a "dual-Meissner" effect. 

Together with the assumption that the effect of the charged off-diagonal degrees 
of freedom are qualitatively irrelevant for the description of the low-energy physics, 
the above scenario is usually called "abelian dominance" . 
Two natural questions arise in the "abelian dominance" scenario: 
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1) Is there a "monopole field operator" which plays the role of an order parame- 
ter for the confinement-deconfinement transition, i.e. with vanishing expectation 
value in the deconfined phase and non-vanishing expectation value in the confined 
phase, so that the transition can be interpreted as due to a "monopole condensa- 
tion"? 

2) does the choice of the field X influence the behaviour of the expectation value 
of the monopole field operator? 

In [6,7] an attempt has been made to give a positive answer to the first question 
on the basis of a circle of ideas which can be presented as follows. 

In electrodynamics with charged scalar field (p one can construct gauge-invariant 
charged field operators and their correlation functions adapting the Dirac recipe 
[8], dressing the local non-gauge invariant field 4>(x) with a cloud of soft photons, 
represented by multiplication by a phase factor with argument given by the gauge 
field A weighted by a classical Coulomb field E(x). 

In abelian gauge theories there is a natural notion of duality exchanging the 
role of charges and monopoles. One can obtain monopole correlation functions 
from gauge-invariant charged correlation functions by a duality transformation; 
in particular this applies to monopoles in U(l) lattice gauge theory. 

In [6,7] it was suggested to apply this construction to the SU (2) gauge theory in 
the abelian projection gauge and it was shown [7] how to render this construction 
gauge-invariant . 

Starting from the expectation value of the monopole operator constructed in 
this manner, Montecarlo simulations show that 

1) this monopole field operator indeed behaves as a good order parameter signaling 
the confinement - deconfinement transition; 

2) the (physical) temperature of the transition is independent of the choice of the 
scalar field X [9]. 

This approach however presents a foundational problem in spite of its great 
numerical success: it is inconsistent with Dirac's quantization condition of fluxes. 
This, in turn, implies an inconsistency of the treatment of small scales. 

To understand the origin of the Dirac quantization condition, let us consider 
the electrodynamics of electrically and magnetically charged point-like particles. 
The equations of motion proposed by Dirac read: 

d»F„ v (x) = q e f v {x\ e^d v F pa {x) = q m j™(x) (1.1) 

where q e jl and q m j u are the electric and magnetic currents generated by the 
particles, q e and q m being their electric and magnetic charges. 
Since equation (1.1) implies current conservation 
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Pjl = o = d»c, 



Poincare's lemma ensures the existence of antisymmetric tensor fields n 
such that 



e 
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However these fields are determined only modulo the transformation 



n 



# 



nt + 8 [p \* 



pa '"pa 

where # = e, m and Xf is a vector field. 

Schwinger proposed an action leading to the Dirac equations (1.1): 



(1.2) 



(1.3) 



\{d { ,A v] + qm n™) 2 {x)+q e {A^){x) 



d 4 



where if the support of the current j™ lies on a curve 7, then n^, should be 
taken as a surface current whose support is a surface E, with boundary 7. The 
current n m describes the Dirac strings attached to the magnetic monopoles whose 
worldlines are given by 7. 

The quantum theory corresponding to the classical action S(A,j e ,n m ) is well 
defined, provided the partition function 



/ 



VAe iS(A,j e ,n m ) 



is independent of the choice of the surface E, with n m satisfying (1.2), i.e. invariant 
under the gauge transformations (1.3). 

If one chooses n 7 ^ and n'™ corresponding to the surfaces E and E', then the 
parameter \™ appearing in (1.3) is dual to a volume-current with support in the 
volume V whose boundary is given by the closed surface difference of E' and E. 
The consistency condition turns out to be given by 



qmqej \™j eil d 4 xe2iTZ. 



(1.4) 



Since are line currents with integral coefficients one can recognize in (1.4) the 
Dirac quantization condition: 



q e q m e 2ttZ. 



(1.5) 



From equation (1.4) it follows as consistency requirement that the integral over 
an arbitrary volume of any electric current appearing in the partition function, 
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multiplied by ^ must be an integer. One can easily prove that this condition 
extends to all physical correlation functions. 

However in the construction of the gauge invariant correlation functions for 
charged fields following the Dirac recipe, one introduces the electric smooth current 
E, whose integral over a generic volume V is a real number, thus violating Dirac 
quantization condition. One may say that in the presence of the electric current E 
the position of the Dirac strings (of monopoles) crossing the support of E become 
visible, thus introducing a physical inconsistency. 

The monopole correlation functions may be obtained from the charged corre- 
lation functions by duality, therefore in the presence of dynamical charges their 
construction following Dirac's recipe encounter the same inconsistency. This ex- 
tends also to the constructions of monopoles for SU(2) presented in [6,7], since it 
involves charged degrees of freedom corresponding to the off-diagonal components 
of the Yang-Mills gauge field. 

Nevertheless, since in U(l) gauge theories the constraint of Dirac quantization 
is expected to become irrelevant at large distances, i.e. in the scaling limit, one 
might still imagine that the low-energy physics, (e.g. the deconfinement transition 
temperature), should not be affected by the above problems. 

In this paper we propose a new order parameter for the confinement - deconfine- 
ment transition, defined in terms of the expectation value of a (regular) monopole 
field operator, which avoids the inconsistency discussed above. 
Our construction is a variant of that proposed for the 't Hooft-Polyakov monopoles 
in the Georgi-Glashow model [10] and, as a basic difference with respect to pre- 
vious constructions, the correlation functions of monopoles are obtained in terms 
of sheets of center vortices. 

The new order parameter exhibits the following features: 

1) it respects the Dirac quantization condition 

2) it is naturally independent of the choice of a U(l) subgroup of SU{2), needed 
in the abelian projection 

3) if one chooses an abelian projection gauge we argue that in the scaling limit 
it approaches the order parameter constructed in [7] in correspondence to that 
projection gauge 

4) it creates a bridge between "abelian" and "center dominance" suggesting how 
one can reconcile the two approaches. 

Although our discussion is heuristic and partly conjectural, the overall picture 
that emerges is consistent with known mathematical estimates and numerical lat- 
tice simulations. 

In order to make the paper selfcontained, we start by recalling in section 2 Dirac 
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recipe and its dual, sketching the corresponding construction of charged and 
monopole Green function in the euclidean formalism. 

In section 3 we review the modification of Dirac's ansatz proposed in order 
to satisfy the Dirac quantization condition if dynamical charges and monopoles 
coexist. 

In sect. 4 we sketch how one defines the magnetic charge in SU(2) theory. 

In sect. 5 we review, in the euclidean formalism, earlier constructions of monopole 

field in SU(2), based on the dual of Dirac's recipe. 

In sect. 6 we define our new monopole field and the corresponding order parameter 
and discuss the link to previous constructions. 

In sect. 7 we outline the connection between "abelian" and "center dominance" 
suggested by our order parameter. 

In order to simplify our formulas in the rest of the paper we use the language of 
forms and currents, both in the continuum and on the lattice. The basic definitions 
of this formalism are reviewed in an appendix. 

2. Dirac's ansatz and its dual 

We start by discussing a simple model, scalar QED, where only electric dy- 
namical charges appear, but no magnetic monopoles. We show how to construct 
gauge-invariant, charged field operators following Dirac's ansatz. 

Let (j) be a massive scalar field with charge e coupled to an abelian gauge field, 
A^, described in terms of a real 1-form A, with classical action 



2 _ 



/rl 1 m l - 

-(«VM) 2 0r) + - \d, - ieA^\\x) + y #W 

or, in the notation of differential forms, 

S(A,d>) = i||cL4|| 2 + I||(d- leA)^ + ^||</>|| 2 . (2.1) 

Dirac's ansatz can be formulated as follows: Let <fi and A denote the quantum 
field operators corresponding to the classical fields <p an d A and let E^dx^ = E x 
denote the 1-form corresponding to the classical electromagnetic field generated 
by a pointlike unit charge located at x in R 3 . Then Dirac's charged field is defined, 
heuristically, by the formula: 

j>(E 3 ) = 4>(x)e^ Aa * E \ (2.2) 

This construction has been rendered rigorous (in the presence of an ultraviolet 
regulation) in [11], within the indefinite metric approach (Gupta-Bleuler gauge). 
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There is however an alternative route to construct charged fields. One starts 
from euclidean Green functions and then invokes the Osterwalder-Schrader recon- 
struction theorem [12]. In approximate terms it works as follows: One constructs 
gauge-invariant euclidean correlation functions for charged fields obeying the O.S. 
axioms (essentially translation invariance, reflection (O.S.) positivity and cluster- 
ing) from which one can reconstruct a Hilbert space of physical states containing 
the vacuum vector Q, a unitary representation of space-time translations, whose 
generators satisfy the spectral condition, and which leaves O invariant, and quan- 
tum field operators. 

There is also a version of the reconstruction theorem that applies to lattice 
theories. 

A euclidean version of Dirac's ansatz is then obtained by replacing the quantum 
field (2.2) by a euclidean field 

(j){E x ) = <f)(x)e J* 4 AA * EX = <j>{x)S A ^ (2.3) 
where E x is the 1-current in R 4 given by 

E x (y) = E s (y)8(y°-x ). (2-4) 

With a lattice regularisation euclidean correlation functions of these fields, 
(Hi (p(Ej ) n, 4>(Ej))i have been proved to satisfy the O.S. axioms [13]. Here 
(■) denotes the euclidean expectation value corresponding to the action (2.1) and 
4>{E X ) is the complex conjugate of (p(E x ). 

It is useful for later purposes to notice that one can obtain a representation of 
the correlation functions of 4>(E) as partition functions of a gas of closed electric 
currents coupled to A, by integrating out </>. 

As an example consider the two-point function 

(cf>(E x )$(Ey)) = i J DAe-^ dA ^det(-A eA + m 2 ) 

e ie{A^- E y )){(j){x) ^ y)){A) (2 _ 5) 

where (-)(A) denotes the (normalised) expectation value corresponding to the ac- 
tion of 4> coupled to A, viewed as an "external" field, and A e A is the covariant 
Laplacian. 

Using a euclidean version of Feynman's path-integral formula for the quantum- 
mechanical time evolution kernel (e sAeA )(x, ?/), one obtains formally 

poo 

(<t>{x)4>{y){A)) = {-A eA + m 2 )-\x 1 y) = / dse~ sm2 (e sA ^)(x, y) = 

Jo 
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J dse~ sm2 J Vq(t)e Jo 



\q\t) + eA^q{t))e{t) 



dt 



q(0)=x 

q(s)=y 

If with every trajectory {(/ M (t), t G [0, s]}, s 6 R + , we associate the 1-current 

3x V (y) = J dtq^(t)S(q(t) - y)dy fl , 
(*jxy is Poincare dual of the trajectory, see appendix), then, one finds 

(cf>(x)^y))(A) = [ v^ 3xy y< A ^\ 



for a suitable measure Vfx(j xy ). 

One can also express the determinant in (2.5) in terms of a sum over closed 
current networks j, so that, for a suitable measure T>n(j) on the current networks 
j, one obtains 

(<f>(E x )t(E y )) = j VAe-^ 2 j VvmK3xy)e ieiMEX - EV+j * y+i)) 

[JvAe-^ dA ^ 2 J P^a)e ie(AJ) ] _1 - (2-6) 

We notice that, while the current networks j in the partition function (the denom- 
inator of (2.6)) are all integer- valued, the currents appearing in the numerator of 
(2.6) also involve a real-valued term (E x — E y ). 

We now sketch how one can obtain monopole correlation functions in a dual theory. 
For [y(l)-gauge theories in d=A dimensions, S-duality is a transformation mapping 
the correlation functions of the original gauge theory onto those of a dual gauge 
theory, exchanging the role of charges and monopoles. The underlying idea can 
be presented as follows. 

Let S(dA) denote a gauge-invariant action for the gauge field A, written in terms 
of its curvature 2-form dA. Introducing an auxiliary 2-form F we can write the 
partition function as 

Z = J VAe~ s ( dA) = J VFe- s ^5(dF) 7 (2.7) 
since the solution of the constraint 

dF = (2.8) 
is given by F = dA and the Jacobian is field-independent. 
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We now express the constraint in (2.8) by a Fourier representation of the 5 func- 
tional: 



S(dF) = JvAe l I FAdA (2.9) 

where A is a new gauge field, the "dual of A" . We define the dual action S(dA) 
through the functional integral Fourier transform 



e 



-S(dA) = J T)Fe -S(F) e ijFAdA^ (21()) 

Plugging (2.9) and (2.10) in to (2.7) we obtain 

Z = J VFe~ s ^ J VAe l I FAdA = J VAe~ § ^ 

where the last term gives the partition function of the dual theory. 

The same procedure proves that duality exchanges the Wilson loop, which can be 

related to worldlines of a charged particle-antiparticle pair, with the Wegner - 

't Hooft disorder operator [14], related to worldlines of a monopole-antimonopole 

pair. 

Let E be a 2-dimensional surface; the Poincare-dual current is also denoted by E. 
The Wilson loop W a (T,), a G R, is defined in terms of the 2-form F appearing in 
(2.7) as 

-ijt ice I dA ice I F ice \ FAT, /o i i \ 

Iy a (E) = e J s = e J £ = e J . (2.11) 
In the same model the Wegner-'t Hooft disorder field is given by 

D a (Z) = e-^ F+ ^- s ^l (2.12) 
The duality transformation acts on such fields as follows: 

(W a (E)) = J VFe~ s ^ J VAe l I FAdA e ia I FAJ: = J V A e -^ dA+a ^ 

= J VFe~^ F+a ^ J VAe l f FAdA = (A*(E))~ (2.13) 

where (-)~ denotes the expectation value in the dual theory. 

To apply duality to scalar QED, we notice that the representation of the partition 

function in terms of the current networks appearing in (2.6) can be viewed as the 

expectation value of a weighted sum of Wilson loops, e Js u> , with a weighting 
measure T>/i(j) in a gauge theory with action 
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S(dA) = ±\\dA\\ 2 , 

if we associate to every current configuration j a 2-current E(j) satisfying 



dm=*y (2.i4) 

As a result the partition function of the dual theory can be written as 

Z = J Vfi(j) J T>Ae-^ dA+e ^ 2 . (2.15) 

Obviously it corresponds to a Maxwell theory with gauge potential A coupled to 
monopoles, whose worldlines are described by j; E(j) can then be identified as the 
surface spanned by the Dirac strings of the monopoles. 

The partition function is independent of the choice of the Dirac strings, since a 
different choice E'(j) also satisfying 



dE'CJ) =* J, 

differs from E(j) by an exact 2-form dV which can be absorbed by a change of 
variables A — > A + V; (Poincare's lemma). 
By performing the shift 



A^ A + e5A _1 E(j) 

and using the Hodge decomposition for E(j) (see equation (A. 4) in the appendix) 
one can alternatively rewrite the partition function (2.15) as 

Z = J j VAe-^ dA+eSA ~ 1 ^ 2 . (2.16) 

The term 5A _1 *j can be interpreted as the magnetic field generated by the mag- 
netic current networks j. One can obtain the two-point monopole correlation 
function, (m(B x )m(B y ))~ , from (2.6) by applying the duality transformation. 
Setting B =* E, and e = g (the magnetic charge in the dual theory), one finds 

(m(B x )m(B y )y = ((p(E x )(p(E y )) = 

i j ^(j)^^)^-^!^^" 1 ^^^-^-^)!! 2 . (2.17) 

The magnetic 3-current B x (y) can be related to the magnetic field strength (2- 
form) B x of a classical monopole located at x in R 3 by 
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B*(y)=B s (y)A5(x -y°)dy . (2.18) 

These ideas can be applied to general models of gauge theories with only dynamical 
charges or monopoles, and they can be rendered mathematically rigorous for lattice 
theories. 

As an example of a lattice theory with monopoles, one may consider the U(l) 
gauge theory in the Villain or Wilson formulation. 

The basic field is a U(l) lattice gauge field 9, i.e. a t/(l)-valued 1-form, and in 
the Villain formulation one must introduce a 27rZ-valued two-form field n. The 
actions are given by 

S v (6, n) =[3\\d6 + n\\ 2 

s w (e) = pJ2( 1 - cos ( de )p) ( 2 - 19 ) 

v 

where the subscript V stands for "Villain" and W for "Wilson". The monopole 
two-point functions are given by 

Gv ( x ,y) = ±. I ft de<xy> J2e-^ de+n+ ^ SA ~ 1 ^ +B ^ B y^ 

V J <xy> n 

G w (x f y) = — [ TT de <xy> e-^- cos (^+^A-( W+B .-B,) p )) j (220) 

Zw J <xi> 

where Z denotes the partition function of the model, the summation in the Villain 
model is over the 27rZ-valued n-configurations and w is a 3-form Poincare dual 
of a path joining {x} to {y}. 

The Green function for monopoles in the Villain model can be recast in a form 
similar to that appearing in eq. (2.17) by defining a real-valued 1-form A and a 
2/rZ-valued 3-form to, with 

A = 9 + <5A -1 n, to = dn. 

Then 

G v (x,y) = -^r- Yl J II dA <*y> 

^ m:dm=0 <xy> 
e -/3||cL4+<5A- 1 (m+27r(B !!; +B J ,+^))|| 2 (2 21) 

Gauge-fixing for A (or quotienting w.r.t. gauge transformations) is understood in 
the A-measure in (2.21). 
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Finally, by denning , e.g. in the Wilson formulation, a disorder field 

D UJ (B X ,B y ) = e -lSw(dO+2nSA- 1 (B*+u 1 -By))-S w (d9)}^ (2.22) 

one can express the monopole two-point function as an expectation value of the 
disorder field: 



G w (z,y) = (D u (B x ,B»)). 

It has been rigorously shown in [13,15] (see also [16]) that this disorder field is a 
good order parameter for the confinement-deconfinement transition in d = 4 U(l) 
gauge theories, G(x, y) approaches a finite value in the confining phase and van- 
ishes in the deconfined phase, as \x — y\ — > oo. 

One can the interpret the non-vanishing asymptotic value of G(x, y) as a signal 
of monopole-condensation. 

3. A modified Dirac ansatz consistent with Dirac's quantization condi- 
tion 

In this section we discuss the modification of the construction of the previous 
section needed when dynamical charges and monopoles coexist. We start by show- 
ing that, as it stands, the above construction becomes inconsistent in this enlarged 
setting. 

Consider, for example, a model of "compact" scalar QED with a scalar field <p 
of electric charge e and monopoles of magnetic charge g, whose partition function 
is given by 

Z = Jvj2(}) J VAe-^ dA +&0)\\ a Jv<pV^e-I^ { - AsA+rn2 ^, (3.1) 

where E(j) are the Z-valued Dirac strings of the monopoles, and the measure 
T>jl(j) is derived from the action of a matter field coupled to A, e.g. a complex 
scalar field of mass rh, through an equation like 

1 ^ rjT { A+ 2)i = / m*)e^ (3.2) 

In order for Z to be physically well defined, it should be independent of the 
choice of the Dirac strings satisfying = j. 

It is convenient to rewrite the integral over and in terms of Z-valued electric 
current networks 1, as 



12 



Then the situation analysed in the introduction emerges, and the consistency 
condition that guarantees that Dirac strings are invisible, in the sense that E(j) — > 
E'(j) is a symmetry of the theory, becomes invariance of Z under the shift A — > 
A + gV(j), for V(j) a Z-valued 1-current satisfying rfF(j) = — E(j), i.e. 

e ie S (va),i) = L (3.4) 

The condition (3.4) is the Dirac quantization condition at the level of currents, 
and it is satisfied, provided the Dirac quantization condition for charges 



eg e 2tvZ (3.5) 

holds. 

If one tries to construct the 2-point function of the charged field according to 
Dirac's ansatz, as in the previous section, one meets an inconsistency, since, in 
contrast to (3.4), 

e ieg(VQ),(l+l xy +E x -Ey)) ^ ^ ( 3 _g) 

even if (3.5) holds. 

The origin of this problem is a violation of Dirac's quantization condition at the 
level of currents, due to the introduction of the real-valued current E. 

Let us first consider this problem for minimal chargesieg = 2n. One might 
envisage avoiding the difficulty encountered above by replacing the electric current 
E x by an electric "Mandelstam string" [17] 7 1 carrying a unit flux along a path 
starting from x and reaching infinity, at fixed time. 

Such a current would still satisfy d^ x = S x , and it does not violate Dirac's quan- 
tization condition. However, the current 7 X does not decay at infinity (in contrast 
to E x ) and, as a consequence, infrared divergences appear in the construction of 
charged correlation functions based on this ansatz. 

In fact, lattice calculations [13] suggest that an abelian gauge theory with massive 
monopoles and charges scales to a gaussian gauge theory at large distances, in the 
Coulomb phase. 

Hence, to every Mandelstam string is associated an infinite positive self-energy, 
~ (7 X , A _1 7 x ), and the interaction energy between two Mandelstam strings 7 X , *y y 
of opposite charge is infinite and negative, ~ — (7^, A~ 1 7 J/ ), because the strings 
have infinite length. 
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Even if the selfenergies are subtructed off, via a multiplicative renormalisation, 
the interaction between the two strings cannot be removed without violating reflec- 
tion positivity, because it depends on the distance between x and y. A violation 
of reflection positivity would, however, render impossible the reconstruction of 
charged quantum fields. 

A possible way to circumvent this infrared divergence was suggested in [10] : one 
has to replace a fixed Mandelstam string by a sum over fluctuating Mandelstam 
strings weighted by a measure Vv(^ x ) supported on strings fluctuating so strongly 
that their interaction energy remains finite even in the limit of infinite length. 
The strings appearing in the construction of the correlation functions of charged 
fields should then converge to a common point at infinity. 

As will be reviewed later, it has been shown in [10] that there exists a lattice 
regularised complex measure T)ve(i x ) on Z-valued 1-currents satisfying 



h x = ^ (3-7) 



such that 

i) the correlation functions for the euclidean fields 



4>(y) J Vv E {l v )e- ie{A ^ V) (3.8) 

satisfy a lattice version of the O.S. axioms; and 
ii) in the scaling limit, 

J Vv E {i x )e ie ^ x > A) ~ e< E ^ A \ (3.9) 

up to a multiplicative renormalisation, where E is the electric "Coulomb" field. 
Thus on large scales, the sum over fluctuating Mandelstam strings reproduces the 
behaviour of phase factor appearing in the Dirac ansatz. (This has been verified in 
[10], in a gaussian approximation). Equation (3.9) suggests that, at large scales, 
the measure Vve(7 x ) mimics an approximate 5— function peaked around E x . 

Let us suppose that an appropriate variant of the O.S. axioms is satisfied by 
expectation values of the euclidean fields (3.8), (as, follows formally, from their 
definition). Then from their correlation functions one can reconstruct quantum 
field operators 
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See [10,13] for details. 

If we consider compact scalar QED and the Dirac quantization condition (3.5) is 
satisfied in the more general form g = ^-q, q £ Z (q ^ 0), then we can repeat 
the above construction of charged quantum fields, replacing the Z-valued currents 
7 X with Z/q-valued currents 7 X (satisfying (3.7)). 

Let us describe how to construct a monopole-monopole Green function. 
We associate a Z/q- valued 2-current £(7 X — ^ y + j xy ) to an integral 1-current 
j xy satisfying 5j xy = 5 y — 5 X and a pair of Z/g-valued currents 7 X and 7 y ( with 
5"f z = S z , z = x, y), such that 

*dE( 7 * _ 7 v + jxy ) = 7 * - 7 » + Jxy . (3.10) 
We then define a disorder field by setting 



£>(£(7 X - <y y +j xy )) = e -^{\\ dA +9{ni x -i v +j*y)+nm\ 2 -\\dA+ni)\\ 2 } , ( 3-11 ) 

One can easily verify that (_D(E(7 X + 7 y + j xy ))) depends only on 7 X + 7 y + j xy 
and not on a specific choice of E, because different choices can be mapped onto 
one another by a shift of A. 

The monopole 2-point function is given by 



/ 



VMl X WMl V )mj X y)(D(m X +1 V + jxy))), (3.12) 

where T>jj(j xy ) is the measure defined through the equation 



— — - = / V^( 3xy )e 



for arbitrary A. 

Denoting by (-)~ the expectation value in the dual model, one can easily verify 
that (3.12) equals 

^(E y )kE x )r- 

From correlation functions of disorder fields such as (3.11) one can obtain, 
via O.S. reconstruction, monopole field operators, m(I? x ), where B x =* E x en- 
codes the infrared behaviour of the soft photon cloud accompanying the monopole, 
(which is an "infra-particle"). 

In particular, for y° < < x°, the correlation function (3.12) is equal to 
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(m(B y )n\m(B x )Q). 

The modification of Dirac's recipe and its dual suggested above can be adapted 
to any abelian theory of coexisting charges and monopoles, and it can be made 
precise for lattice theories. The measure Due can then be constructed as follows; 
(we omit some technical details about b.c, see [10 ]): We introduce a £7(l)-valued 
scalar field x °f period 2-7rg defined on a three-dimensional hyperplane A x o, at 
fixed time x°, minimally coupled to the gauge field A. An action is given e.g. by 

S( X ,A) = I £ [l-oos ( d * + -*)<*»> ], (3.13) 

<xy>eA x0 

and we denote the corresponding expectation value by (-) x o(A). 

Then the two-point function, (e lXx e~ lXv ) x o (A), can be expressed in terms of 
Z/g-valued 1-currents ~j xy satisfying 5"f xy = 5 X — 5 y as 

(e^e-^y) x o(A) = J Dv{ lyx )e< A ^\ (3.14) 

where x° = y°. 

If £ is sufficiently large (and if the field strength {(cL4) p } pe A^ does not fluctuate 
much) the system described by x is i n a phase where the symmetry 

X -> X + const (3.15) 

is spontaneously broken [3]. The associated Goldstone boson is a real field, A, 
describing a spin wave. It corresponds to the decompactification of the range of 
X- Deviations from the theory described by non-interacting spin waves are due 
to vortices. They are believed to be irrelevant at large scales (in the sense of the 
renormalization group), provided the symmetry (3.15) is spontaneously broken. 
We propose to evaluate (3.14) at large scales, neglecting vortices. The result is 

Z ■ ™e«. e -«.e-^»— II 
= e -|((^-' 5 y). A ;o 1 (^- ,5 s')) e *( dA ^ 1 (^-^).^) 5 (3.16) 

where £ ren denotes a renormalised coupling, and A x o is the 3D laplacian on A x o . 
One can easily verify that 

E*-E» = dA£(6 x -6 y ). (3.17) 
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Sending y to oo and eliminating the first factor in (3.16) by a multiplicative renor- 
malization, one obtains 



This suggests that the desired measure is of the form 



Vv E (l x ) ~ lim c y Vv{^ X y), 



where c y denotes a suitable multiplicative renormalisation constant, that we expect 
to be finite, on the basis of a gaussian computation performed in [10]. 

As example of lattice theories with dynamical charges and monopole one may 
consider models where a {/(l)-gauge field 9 is coupled to a £7(l)-valued scalar 
matter field <fi of charge q = 1,2,...; (here we set the elementary charge equal to 
unity) . 

In Wilson's formulation, the action is given by 

S = - cos(d0) p ) +k^(1- cos(g0 + d0)<y>). (3.17) 

P <ij> 

According to the recipe explained above, the monopole two-point Green function 
(for monopoles of unit charge) is given by 

\ \ n d»<ii>n^* / vmi x )vmi v ) 



■^E P ( 1 - cos ( de+27rE ^ :!; -^ !/ +^^)p) P - K E< iJ >( 1 - cos ('? e +^)<^>) 



(3.18) 



where £(7 X — 7 y + ? xy ) is a Z/g-valued 2-form satisfying (3.10), and j xy is an 
integral 1-current with support on a path connecting x to y and Jjxy = 5 X — S y . 

The Green function (3.18) is independent of the choice of j xy since a change 
in j xy can be compensated by a shift of 9. Hence, in contrast to the continuum 
models previously discussed, no integration over j xy is needed. 

4. The magnetic charge in 577(2) gauge theory 

We turn to the SU(2) Yang-Mills theory, identifying its magnetic monopoles 
and defining their magnetic charge. 

We first discuss these matters in the continuum, where we denote by A the 577(2)- 
connection form. 

Let us consider a cube c in R 3 ; its boundary, dc is homeomorphic to a 2-sphere. 
The restriction of A to dc can be viewed as a connection form A = A\g c of an 
50(3)-bundle over dc. 
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S'0(3)-bundles over S 2 ~ dc are classified by iri(SO(3)) ~ Z 2 . The following 
relation holds between A and the integer class n mod 2 e Z/2Z ~ 7r 1 (S'0(3)) 
classifying the corresponding SO (3) bundle: 

where A( p ) is the connection 1-form of an S't/(2)-bundle over the face ("plaque- 
tte") p of the cube c, obtained by lifting the 5'0(3)-bundle over p. 

If e lnn 7^ 1 one cannot extend the 5'0(3)-bundle to the interior of the cube c; 
this signals the presence of an odd number of Z 2 -monopoles of SO (3) inside c and 
we identify e mn as the Z 2 -charge contained in c. 

An abelian projection gauge, defined as in the introduction, selects a £7(1) resid- 
ual gauge group. By projecting A to the Cartan subalgebra of su(2) corresponding 
to the residual gauge group one obtains a £7(1) connection a on dc. 
£/"(l)-bundles over S 2 ~ dc are classified by ni(U(l)) ~ Z. The relation between 
the integer n classifying the bundle and the connection a is given by 




where F(a) is the curvature of a. 

We identify the magnetic charge contained inside c in the abelian projection 
gauge with the integer n, corresponding to the first Chern number of F(a). 

If a is derived from A by abelian projection, the integer n appearing in (4.2) is 
the same appearing (4.1). 

According to a general theorem [18], the classical monopoles associated to a 
£7(1) subgroup SO(3) are regular if they carry even magnetic charge. For these 
monopoles there are gauge choices for which no singularity of Yang-Mills curvature 
occurs at the monopole position. Monopoles with odd magnetic charge are singu- 
lar, i.e. in every gauge there is a singularity of the Yang-Mills curvature where 
the monopole is located. The position of the singularity is then independent of the 
choice of the £7(1) subgroup, or equivalently of the abelian gauge projection, used 
to define the magnetic charge and it identifies the position of a Z 2 -monopole. 

The definition of magnetic charge and Z 2 - charge can be easily adapted to the 
lattice formulation as follows [19]. 

Let g denote the ££7(2) lattice gauge field. The Yang-Mills action on the lattice 
is given by 
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Sym(9) = (3J2( 1 -x(98p)) (4-3) 

v 

where x is the character of the fundamental representation. 

Denote by X(g) the scalar field with values in su{2) transforming under the 
adjoint representation, identifying the abelian projection gauge. Let W be the 
SU(2) gauge transformation such that for every site i on the lattice one has 

W i X i (g)W t + = \X i (g)\<r 3 , (4.4) 
were \Xi\ = ^TTXfXf and define 

9<ij> = W l9<lJ> W+. (4.5) 

g denotes the S , L r (2)-gauge field in the abelian projection gauge defined by X. 
g<ij> can be decomposed as a product of two matrices C<ij> and u < ij > (9), where 



(1 — |c<ij>| 2 ) 2 —C<ij> \ „, _ i\ 



c <->-{ c<ij> (i-ic< y> nU' tt <->- e2 ' (4 - 6) 

with c<ij> G C, c<ij> denoting its complex conjugate, and \0<ij> = arg(^<ij>)n- 
Hence 6> is a t/" (1) gauge field with range (— 2/T,27r) and c is a charged field of 
charge 1. 

A similar decomposition holds for the coset variable given on a link < ij > by 
(7<ij>r = U < ij > , where Y is the centre of SU(2), isomorphic to Z2. U<ij-> can be 
viewed as an 5'0(3)-gauge field and in the decomposition (4.6) 9 is now a U(l) 
field with range (— tt, it). 

We define the magnetic charge in a lattice cube c by 



pEoc 

where (d9) p is restricted to the range (—tt, tt). 
The Z2-charge in c is defined by 

e i*z c (g) = e l E pe ac ar S^^). (4. 8 ) 

A plaquette p where 
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can be identified as the location of a Dirac string of a Z 2 - monopole intersecting 
the plane containing p. 

In the lattice formulation the relation discussed in the continuum between Z 2 - 
and magnetic charge becomes 

e iirz c (g) _ e iirm c {9) (4-10) 

5. The U(l)— monopole order parameter for SU(2) Yang— Mills theory 

In this section we recall some attempts made to derive an order parameter in 
SU(2) Yang-Mills theory on the lattice, using a construction directly inspired by 
Dirac ansatz. 

Integrating out the charged field c defined in (4.6) one can view the SU(2) 
gauge theory in an abelian projection gauge as a U(l) theory with gauge field 
9 and an effective action of the form S e ff(d0), where we used the residual U(l) 
gauge invar iance to deduce the dependence on the curvature d9. 

It is then natural, following the dual of Dirac ansatz, to try to define the two- 
point Green function of the magnetic monopole of the abelian projection by per- 
forming the shift 

d9^d9 + 47T5A" 1 (B x -B y + u) (5.1) 

where w is a 3-form Poincare dual of a path connecting x to y. 

This is basically the attempt made in [6] in the Maximal Abelian Gauge [19] 
and in [10] in the Spatial Maximal Abelian Gauge. 

However, a closer look shows some inconsistency in this approach. In fact, it is 
easy to prove that the effective action is given by an expression of the form 

J2c c e qa3 ^<^>^ 9<lj> (5.2) 
c 

where £ is a loop and Cjr a complex coefficient independent of 9. 

Choosing, for each £, a surface E(£) such that <9E(£) = £, one can rewrite (5.2) 

as 

c 

i.e. in the form S e ff(d9) and the choice of the surfaces £(£) is obviously irrelevant. 
However, when the shift (5.1) is performed the term 



(5.3) 
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is no more independent of the choice of £(£). This inconsistency is due to the 
violation of Dirac quantization condition at the level of currents related to B x —B y . 

The approach of [7,9] is slightly different and monopole Green functions are 
defined directly by modifying the Yang-Mills action. 

Let us fix a convention assigning to every plaquette p a site j(p) on its boundary 
(see [9] for more details) and let X(g) denote the scalar field defining the abelian 
projection. 

Within our setting the proposal made in [7,9] to obtain the monopole-monopole 
Green function is to replace the plaquette term x(ddp) m the Yang-Mills action 
(4.3) by 

x ( gdpe ^U P )MSA-\B*-By + „ )]p y (5 _ 4) 

Actually, in [7,9] the support of to consists of a sum of the dual of two straight lines, 
uj x and — uj v each at constant time in the 3-direction. With this choice B x + uj x 
can be related, as in (2.18), to the magnetic field of a Dirac monopole with its 
Dirac string along the 3-direction. 

In the abelian projection gauge defined by X, using the decomposition (4.6), 
one can rewrite the argument of x m (5.4) as 



( e ic<ij>-a e i±9 <ij> <T3^ e i<T 3 2TT[5A 1 (B x -B v +uj)] p ^ 

<ij>€dp 



where 



c <ij> — C <ij> ^ C <ij>i C <ij>i ^ a — 1> ^ C — (c , C ). 

We insert the lattice constant e and, since we are really interested in the continuum 
limit, we consider an expansion in e up to 0(e 2 ): One finds 

(5.5) ~ fl (e iS<i i>- 3e ) e^s^ide+SA-^B^-By+u;)]^^ + o( e 2 )). 

<ij>Edp 

Therefore to 0(e 2 ) this recipe seems to give a consistent prescription; however 
problems arise at order e 2 because in gg p there appear terms depending on 9 which 
cannot be rewritten in terms of (d0) p : 

0(e 2 ) = e 2 {[(c$) p a 3 , 7(c) • «?] + (9a 3 A 7(c) • a) p } + 0(e 3 ), (5.6) 

where A is the wedge product on the lattice (see e.g. [10]) and 7 = (7 1 ,7 2 ) are 
functions of c. 
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As a consequence of the last term in (5.6), a change in the choice of the "Dirac 
string" u cannot be eliminated by a shift in 9, and again this can be traced back 
to a violation of Dirac quantization condition at the level of currents. 

The consistency to order e 2 of the recipe in [7,9] suggests however that this 
disorder parameter (and the one of [6]) could be meaningful at large scales, as 
supported by numerical results for the (physical) temperature of the deconfinement 
transition. 

6. Green functions for regular monopoles in SU(2) lattice gauge theory 

Combining the ideas of the last three sections, one is lead to propose a definition 
of Green functions for regular monopoles of any abelian projection, as follows. 

Imagine that a regular charge-2 monopole in an abelian projection gauge is 
created at a site x and annihilated at a site y. We propose to construct the corre- 
sponding two-point function summing over a pair of fluctuating strings carrying 
magnetic flux 1 with end points at the position of the monopole. 

As discussed in section 4, strings of odd magnetic flux in every abelian projection 
can be identified as Dirac strings of Z2-monopoles, and hence are independent of 
the choice of the projection. 

In turn, Z2-monopoles can be introduced by means of a 't Hooft disorder fields 
which, in SU{2) lattice gauge theory, is defined as follows. 

Let E be the Poincare dual of a surface bounded by a loop C. Then the corre- 
sponding disorder field is defined by 

£>(£) = e - [SYM({99pe ^ Sp ° 3}) - SYM{{9ap})] (6.1) 

where Sym is the action (4.3). The expectation value of D(E) depends only on C 
and describes the worldlines of a Z2 monopole-antimonopole pair. 

Since the regular monopoles have even magnetic charge and the charged field, c, 
of the abelian projection has integer electric charge, one can adapt to the present 
setting the construction presented in section 3 for q = 2. 

Let j xy be an integer 1-current satisfying Sj xy = 5 y — 5 X with support on a path 
connecting x to y, and let Vi/e('J x ), T>i/E('j y ) denote normalized, signed measures 
over Z/2-valued currents 7 x ,7 y , constrained by 5"f z = Sz, z = x,y as defined in 
section 3. 

To a configuration {7 X , 7 y , j xy } we associate a 2-current £(7 X — ^ y + j xy ) sat- 
isfying 
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where the 2 in the.h.s. appears for agreement with definition (6.1). E is then 
the Poincare dual of a two-sheet surface, with boundary given by the support of 
7 X — 7 y and the two sheets joining each other along the support of j xy . 
We define the monopole two-point correlation function by 

G(x,y) = J VMl x )Vv E {l y ){D{ni x -l v + 3*y)))- (6-2) 

We claim that long-range order in this correlation functions characterizes the 
confinement phase. 



From the n-point monopole correlation functions defined by generalizing eq. (6.2) 
one can (at least formally) reconstruct, via O.S. lattice reconstruction, a monopole 
field operator rh{B x ),B =* E. The long-range order for G{x, y) then corresponds 
to a non-vanishing vacuum expectation value of rh. 

Since equation (6.2) does not involve any abelian projection, it follows that, 
while the definition of the trajectory of a regular monopole requires choosing an 
abelian projection, the locations of creation and annihilation of a monopole are 
independent of that choice, hence intrinsic to the SU (2) theory. 



From numerical simulations one may gain some indirect support for the conjecture 
that m is a good order parameter for the confinement-deconfinement transition 
by noticing that the large distance behaviour of G(x, y) appears to approach that 
of the £7(l)-monopole Green functions of [7]. 

In fact, since the group manifold of 577(2) is isomorphic, via the exponential 
map, to a 3-ball of radius n with boundary points identified, one may replace 
e i7rE p cr 3 ky e l7r (— )j(p) e p j n definition (6.2) of the disorder field for any choice 
of scalar X defining an abelian projection. This is because (-[x[)j(p) defines a unit 
vector in su(2). 

As remarked in section 3, at large scales the measure Vi> E (^ x ) behaves as an 
approximate Dirac measure peaked around the current configuration E x , and this, 
in turn, implies that in (6.2) the configurations of dY>(^ x — 7 y + j xy ) are peaked 
around 2(B X -BV + with B =* E,u =* j xy . 
Therefore, in the scaling limit, one expects that 

J vMi x )vMi y )e~ SYM ^ e " ( * )j(p)Sp >] 

^-Wto/m',.^'- 1 — (6 _ 3) 
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This would reproduce the behaviour of the disorder field defined via eq. (5.4). 
Numerically, its expectation value provides a clear signal of the confinement- 
deconfinement transition [7]. Eq. (6.3) yields an explanation of the numerical 
evidence that, from order parameters corresponding to different abelian projec- 
tions, one obtains the same transition temperature [9]. 

In fact, this transition is governed by the low-energy physics, correctly captured 
by the scaling limit. Hence, assuming that eq. (6.3) holds, the order parameters 
defined via eq. (5.4), unphysically dependent on the choice of the "Dirac string u" 
at small scales, are just the scaling limit of the order parameter defined via (6.2), 
which is manifestly independent of the choice of an abelian projection and of the 
"Dirac string j xy n . 



Expressing {D\Yi(~f x — 7 y + j xy )j) in terms of magnetic currents, through a 
duality transformation, one can exhibit more explicitly the non-vanishing vacuum 
expectation value of m as a dual Higgs mechanism (in the spirit of the "abelian 
dominance" scenario). 

Since, however, £(7 X — 7 y +jxy) can be interpreted as a (double) sheet of center 
vortices, a connection with the "center dominance" scenario emerges, as discussed 
in next section. 

To derive the duality transformation, we express the correlation function (D(E)), 
where E = £(7 X — 7^ + j xy ), in terms of the variable 6 and c defined in eq. (4.6), 
and insert in the integration measure an abelian projection gauge fixing 



First we integrate out c. As a consequence of the residual U(l) gauge invariance 
we can expand 



II ^<y>*<y>e- 5 ^«^ (fl - c)ei * Ep ' 8 > ) n<4q(^ c ) - ( 6 - 4 ) 



S(—(c,e)-a 3 ). 
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as a Fourier series in dd + 2/tE. 

The coefficients are denoted by F(n), where n is a Z/2-valued 2-form. 
We define a 1-form £ by 



5n = £, 



(6.5) 



and decompose the 2-form n as 



n = n[£] +* d£, 



(6.6) 
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where n[£] is a Z/2-valued solution of (6.5) and £ is a Z/2-valued 1-form in the 
dual lattice. Then one obtains 

(0(E)) = ^EE +* rf / II de <tJ> e^e^^\ (6.7) 

[£] <?:<«=0 J <ij> 

where [£] denotes a gauge equivalence class of £, and the equation 

(nfl,d0) = &0) 

has been used. 

Integrating over 9 imposes the constraint t = 0. Hence, in particular, one can 
choose n[£] = 0. Furthermore we can replace £ by a real-valued 1-form A by 
inserting the term J2 p . Sp=0 e l ^ 1T ^ A,p \ where p is a Z-valued 1-form (see e.g. [21]). 

As a result we obtain 

(D^* +f y ))) = \ f d[A] F{*dA)e lA ^-^ v+ ^y +p ' A \ (6.8) 

p:Sp=0 

where d[A] denotes formal integration over gauge equivalence classes of A. 

In (6.8), worldlines of regular monopoles of the abelian projection are described 

by the currents j xy + p, and they exhibit sources at {x} and {y}. 

The representation (6.8) explicitly proves independence of the choice of j xy in the 

construction of Green functions. 

By setting 

S(A) = - In F(*dA) 

one can view the dual model appearing in (6.8) as a Higgs model with gauge action 
S(A) and the correlation function 

can be viewed as the two-point function of the charged field, (4>(E x )(f)(E y )}~ , of 
that dual model, where (-)~ denotes the corresponding expectation value. 

The abelian Higgs model in four dimensions has two phase, the Coulomb and 
the Higgs phase. If the dual model is in the Higgs phase one expects that 
((fi(E x ))~ 7^ 0. (In fact, for the standard gauge action and with the Dirac recipe 
for the charged field, this has been proved in [13], [16].) In the original model this 
non-vanishing expectation value corresponds to (Q\rh(B x )Q) ^ 0, thus describing 
monopole condensation. 
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The suggestion that the dual model is in the Higgs phase comes from numerical 
simulations, as previously discussed. 

The above construction makes the relation between the original SU(2) gauge the- 
ory and a dual abelian Higgs model more precise, as advocated by many authors. 

We end this section with a remark about about monopole Green functions in the 
formal continuum limit. Our discussion of monopole in SU (2)-Yang Mills theory 
was performed in the lattice, because it heavily relies upon the 't Hooft disorder 
field which has no simple continuum analog. 

Presumably one can construct a disorder field with the desired properties in 
the continuum using the loop space formalism developed in [22], but this will be 
discussed elsewhere. 

7. Relation with center dominance 

The representation (6.8) does not exhibit center vortex sheets; they are hidden 
in the definition of F(*dA). To exhibit them explicitly one starts, following [23], 
by replacing the SU (2)-gauge field g with a couple of new variables {U, a} where 
U is the gauge coset variable, introduced in sect. 4, which can be viewed as an 
SO (3) gauge field and a is a 2-form with values in {0, 1} ~ Z2. 

It has been shown in [22] that the two fields U and a are not independent. It is 
easy to show that J2 P edc ar SX(#9p) is orn y a function of the coset field U, which 
we denote by z c (U) and the following constraints holds: 

e iir(da) c _ e iz c (U) (7.1) 

for each cube c. 

Let us discuss the relation between a and the 't Hooft disorder field. 

The plaquette term x((7a p e l7I " <T3Ep ) appearing in (D(E)) is rewritten in the new 

variables, as 

\x\(U dp )e^+^\ (7.2) 

where \x\(U) = \x(g)\- 

Hence the introduction of the disorder field D(E) induces a shift of a p by £ p . 
The constraint (7.1) can be solved by 

e ina p = s i gnx (C/ ap ) Y[ sign X (t/ < y > ) 
<ij>€dp 

= signx(gdp) JJ signx(0<ij>) (7.3) 

<ij>€dp 
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This solution is gauge-dependent but does not involve any choice of abelian pro- 
jection. 

In the center dominance scenario one defines the maximal center gauge in SU (2) 
gauge theory as the gauge which brings the link variables {g<ij>} as close as 
possible to the center, T ~ Z 2 , of SU{2), by maximizing the quantity 

R = T r(g<ii>) 2 - 

<ij> 

In the maximal center gauge, a plaquette p where 

II sign X (<7<y>) = -1 (7.4) 

<ij> 

is the location of a P-vortex. 

It has been rigorously proven in [23] that, for large f3, i.e., close to the continuum 
limit, the set of plaquettes where sign x(ddp) = — 1 is dilute. 

Therefore, since sign x(ddp) is gauge-invariant, for large /?, the identification of 
a P-vortex location as the set of plaquettes where a p ^ in the center projection 
gauge, should be equivalent to the standard definition, equation (7.4), from the 
point of view of discussing the deconfinement transition. 

Numerical simulations shows that, in the confining phase, P-vortex sheets per- 
colate [2]. This suggests that in this phase even the introduction of an additional 
infinite P-vortex sheet E, like the one involved in the construction of the monopole 
Green function, should be a small perturbation and should not lead to a clustering 
behaviour. 

In other words one could interpret a non-vanishing expectation value of the 
monopole operator m in the maximal center gauge as due to a condensation of 
P-vortex sheets, in the spirit of center dominance. 

In the deconfinement phase at positive temperature, numerical simulations shows 
that P-vortex sheets are dilute [2]. Hence, for large /?, it is natural to conjecture 
that the introduction of an infinite P-vortex sheet S leads to clustering and, as a 
consequence, to a vanishing expectation value for rh. 

Finally we remark that an approximate relation between P-vortex sheets and 
regular monopole worldlines can be established following [24]. 
The double-sheet P-vortex structure associated to monopole worldlines appearing 
in [24 ] is a natural counterpart of our construction of monopoles in terms of the 
double-sheet surface S. 
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Appendix 



Forms and currents in the continuum and on the lattice 

In the continuum we consider as " space-time" the euclidean space R d . 

Given an antisymmetric tensor field of rank k on R d , a Ml ... Mfc (x) one defines the 

associated /c-form by setting 

a ik) (x) = ^a^.^ k (x)dx^ A ... A dx»* (A.l) 

where A is the wedge (antisymmetric tensor) product. The space of /c-forms is a 
group A fc (R d ) under pointwise addition. We denote by 

d: A k (K d ) ^A k+1 (K d ) 
the exterior differential defined through 

da(k H x ) = /i.L. ^i.-ftW^ A dx^ A ...dx»\ (A.2) 

by 



* : A fc (R d ) -> A d " fc (R d ) 



the Hodge star defined through 



1 1 



*(a ( %)) = A A ^ d " fc ' (^- 3 ) 

by 5 =* d* (-1)^+!) the codifferential and by 

A = 5d + d5 (AA) 

the Laplacian. 

An inner product between /c-forms is defined by setting 

(a (k \b^ k) ) = J d d xa^^ k (x)b^-^(x) = J a<*> A* &<*> (A5) 
and it satisfies 

(a (fc) , db {k - 1} ) = (Sa {k) , b {k - 1] ). (A6) 

The L 2 norm corresponding to the inner product (A. 5) is denoted by | ||, i.e. 
(a( k \a^) = \\a^\\ 2 . 
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The Poincare lemma states that if da^ = 0, then there exist a^ k ^ such that 
a (k) =da (fc-i). 

A /c-current in R d is a linear functional in the space of d — k forms with compact 
support, continuous in the sense of distributions, i.e. /c-current are /c-forms with 
distribution-valued components [25]. 

In the space of currents there exist a map, Poincare duality, associating to a 
/c-dimensional surface T, k a (d — /c)-current, PD(T, k ), according to 

f a ^ = [ a,™ ARD(£ fc ), (A. 7) 

for any k— form of compact support. The following property holds: 



PD(dY, k ) = dPD{Y, k ) (A8) 
where d denotes the boundary operator. 

A basic consequence of Poincare's duality is that, whenever well defined, 

/ PD(Z k ) A PD(Z d _ k ) 

is an integer counting the intersection with sign of T, k with Tid-k- 
Linear combinations of such /c-currents PD(E k ) with integer coefficients are called 
integral k— currents. Poincare's lemma holds also for currents. 
We now turn to the lattice. 

Our lattice is Zf , 2 , where the subscript 1/2 indicate that the coordinates of the 
sites are half-integer. 

If W is an additive abelian group, one can define /c-forms with values in W as maps, 
a ( fc ) ? from oriented /c-cells, c k , of the lattice to W satisfying a^ k \—c k ) = —a^ k \c k ), 
where — c k denotes the cell obtained from c k reversing the orientation. 
We denote by d the lattice exterior differential: 

da^(c k+1 )= a(fc) ( Cfc ) ( A9 ) 

Ck€dc k + 1 

and by * the Hodge star. Let c* d _ k denote the cell in the dual lattice, Z 4 , dual to 
c k . Then 

*(« (fc) )(^)=« W ^). (A10) 

We also introduce the (lattice) co differential 5 = (-) d ( fc+1 )*cP and Laplacian A = 
d5 + Sd. 
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If W is a Hilbert space with inner product ( , ) one can define a inner product 
among VF-valued /c-forms and by 



( a ( fc ),6( fc )) = ^(a«,6W). 



(All) 



The £ 2 -norm corresponding to this scalar product is denoted by | ||. Equation 
(A6) and the Poincare lemma hold also on the lattice. 

If W is a discrete group we call the VF-valued k— forms also k— currents, in analogy 
with the continuum definition. 

If Efc is a k— dimensional surface in the lattice one defines its Poincare dual as the 
d — k current PD(Ek) in the dual lattice such that 



In the paper we do not make use of the symbol PD, introduced in this appendix 
for sake of clarity, and often identify a /c-surface with its Poincare dual defined 
above. 
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